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Abstract 

We study open bb production at large rapidity at the LHC in an attempt to pin down 
the gluon distribution at very low x. For the LHC energy of 7 TeV, at next-to-leading 
order (NLO), there is a large factorization scale uncertainty. We show that the uncertainty 
can be greatly reduced if events are selected in which the transverse momenta of the two 
.B-mesons balance each other to some accuracy, that is + p 2 rl < ^0- This will fix 
the scale fip ~ ko, and will allow the LHCb experiment, in particular, to study the x- 
behaviour of gluon distribution down to x ~ 10~ 5 , at rather low scales, \i ~ 2 GeV. We 
evaluate the expected cross sections using, for illustrative purposes, various recent sets of 
Parton Distribution Functions. 



1 Introduction 

The data from HERA and the Tevatron do not constrain the behaviour of the low x gluon 
density, g(x,Q 2 ). Indeed, if Q 2 ~ 4 GeV 2 , then already for x < 10~ 3 there is a significant 
difference between the gluon distributions found in the different global PDF analyses. On the 
other hand, this is just the region sampled by the underlying events at the LHC. 

It appears attractive to use the inclusive bb production at the LHC to study the behaviour 
of the gluon distribution in the very low-x region. Indeed, due the rather large mass of the 
b quark, the process may be described in the framework of perturbative QCD. The dominant 
contribution arises from the gg — > bb hard subprocess. Its cross section has the following 
structure 

da/d 3 p = / dxidx 2 g(xi,fi F ) \M(p; i~l f , Hr)\ 2 g(x 2 , Hf) , (1) 



where the gluon densities, g(xi,/j,p), are taken at some factorization scale //p, and the matrix 
element squared, \A4\ 2 , describes the cross section of the elementary gg — > bb subprocess. The 
process samples gluons which carry momenta fraction x, of the initial protons, where 



m hard 



exp(±y). 



(2) 



Xl,2 = 




At the LHC energy of y/s = 7 TeV, and rather large rapidit}o, y ~ 5, of the whole system 
produced in hard subprocess, one can probe the gluon densities with x ~ 10 -5 . For this estimate 
we have taken the mass created in 'hard subprocess' = 10 GeV. Recall that at present 
there are no data in this small x domain and different global parton analysis predict quite 
different gluons, especially close to the input scale for parton evolution, see, for example, [2j[3]. 
It therefore appears that the LHC, and the LHCb experiment in particular, offers a golden 
opportunity to make a precise determination of the gluon in this important low x domain. 
However, first we must face the problem of the choice of factorization and renormalization 
scales. 

A factorization scale \ip is needed to separate the contributions hidden in the incoming 
Parton Distribution Functions (PDFs) from those that included in the hard matrix element 
| _A^f | 2 . Here, the gluon g(x,fip) is the PDF that we are concerned with. Contributions with 
low gluon virtuality q 2 < fMp are included in the PDF, while those with q 2 > \j? F are assigned 
to the matrix element. The second scale, the renormalization scale /ir, in ([T]) is necessary to 
fix the small value of QCD coupling, a s and to justify the perturbative QCD approach. 
In principle, if all contributions (NLO, NNLO, etc.) are included, then calculated cross section 
would not depend on the values chosen for both of the scales and n F - 

2 Problems associated with the choice of scales 

However, one faces difficulties in the description of the new LHC data [1] for bb production. 
We list these below. 

• The NLO QCD prediction strongly depends on the choice of factorization scale, see Fig.[TJ 
For example, the result obtained with the choice /ip = 2m± is more than twice larger 
than that for the case of jip = m±/2, where here m± = \Jp\ + m 2 ). 

• Moreover, at the NLO, we have a sizeable contribution from the 2 — > 3 (gg — » bbg) 
subprocess, where one additional gluon is emitted in the hard collision. This leads to 
a considerable smearing of the x domain where we sample the incoming gluons. The 
smearing is especially strong if we adopt a low factorization scale, because then there is 
a large phase space allowed for gluon emission from the matrix element. Note that the 

1 Rapidities in this range are optimal for the LHCb experiment [T]. 



2 



60 




LHCb > 







2 



3 4 5 



6 



Figure 1: The NLO predictions for the cross section of bb production obtained using the FONLL 
program [4] from MSTW08 |2J (continuous curves) and CT10 |3J (dashed curves) parton sets, at 
the LHC energies of 7 TeV and 14 TeV, as a function of pseudo-rapidity t)b with scale jjp = mi 
and m b =4.7b GeV; compared with LHCb data at 7 TeV [lj. The predictions using MSTW08 
partons are also shown for four choices of factorization scale: [ip = 2m±, m±, m±/2, m±/4. The 
renormalization scale is set to hr = m±. 

probability of emission is enhanced by two large logarithms^: hi(jn\/ /Mp) and ln(l/x). In 
particular, Fig. [2] shows that if we choose a scale \ip = m_i_/4 then the major contribution 
comes from x ~ 10 -2 , and not from x ~ 10 -5 as we had hoped. From this viewpoint it 
would be better to take a large /ip. 

• On the other hand, to differentiate between the low x gluons it would be better to work 
with a relatively low fip, where the difference between the different global PDF analyses 
is larger. At high scales fip, a large fraction of low x gluons comes from a region of much 
largecl x in input distribution, where the input distribution is already well constrained 
by existing data. Therefore, at larger fip, predictions for LHCb bb production, based on 
different PDF sets, become close to each other. 

2 The ln(l/x) enhancement is the main origin of the scale uncertainty observed in the collinear NLO approach 
at very small x. If we were to decrease the factorization scale /i^, then we have to move gluons with p g x ~ (J,f 
from the PDF to the matrix element. The problem is that, at very low x, there may be several gluons emitted 
in the PDF, while only one gluon emission is allowed in the NLO matrix element. This spoils the compensation 
between the variations of \A4\ 2 and the PDF, which should provide (and, indeed, in the larger x region, does 
provide) the stability of the results under scale variations. 

3 Recall that a parton loses x during DGLAP evolution. 



3 




£3 I I I I I I I I 

g 10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 1 
o 



X2 

Figure 2: The distribution of the values of x 2 of the gluons sampled in NLO bb production with 
Tjt, = 5.5 at the LHC energy of 7 TeV, after the cross section has been integrated over x\ and 
pr- The distributions are shown for the MSTW08 [2J and CT10 [3J parton sets for three values of 
the factorization scale, namely /xp = m_|_ (upper pair of curves), m±/2 and mi/4 (lower pair of 
curves at small x). Note that there is a difference between y and 77, namely 77 — y ~ ln(mj_/pr). 
This difference is accounted for in the FONLL calculation [4] used to produce this plot. The 
renormalisation scale is taken to be /Jr = m±. 



• Moreover, recall that since for a large scale _n F ~ 2m± up to the half of the cross section 
originates from rather heavy virtual gluonsj, the original perturbative calculation, which 
assumes that the gluon virtuality q 2 is small in comparison with the quark mass (or mi), 
becomes inconsistent. 

• Finally, at NLO, we also have an unavoidable uncertainty in the prediction of bb pro- 
duction arising from the choice of the renormalization scale, fiR, which we will discuss in 
Section 

However, despite these difficulties, we show that it is possible to use the LHC bb data to make 
a measurement of the shape of the gluon PDF in the interval 10~ 5 < x ^ 10~ 2 . 



3 Fixing \ip by a cut on the vector sum of px' s 

To overcome the problems associated with the choice of renormalization scale, we may ask for 
the measurement of the cross section of bb events in which the two quarks balance each other in 

4 We have seen from Fig. [T] that <j(hf = 2m±) > 2cr(fj,F = m±/2). 
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the transverse momentum plane to some accuracy; that is p 1T ~ ~P2t- I n other words, to seek 
events which satisfy a cut on the vector sum of the transverse momenta of the heavy quarks, 

\p 1T + p 2T \ < k . (3) 

Of course, in this way we will lose some of the cross section, but this should not be a problem 
with the available high LHC luminosity. Another point is that one cannot measure the quark 
momentum directly. However the momentum of the 5-meson can be measured^], and due to 
the strong leading effect^ in B- meson production, the event selection, proposed in (j3J), can be 
performed with sufficient accuracy for those events with reasonably small transverse momenta 
of the 5-mesons, say, \Ptb\ < 5 GeV. Moreover in order to better constrain the x values of 
the gluons in our selected events, we may put an additional cut on the pseudo-rapidities of the 
S-mesons, say, 

— ffel < 1 or V ~ 0.5 < rji , 772 < r] + 0.5. (4) 

Here, for illustration, we use the latter cut and present results as a function of 77. 

With the above kinematics it is natural to choose \xp = k . At first sight, in this way we 
appear to have excluded any gluon emission due to the NLO matrix element; a gluon with 
a transverse momentum, p g T, less than [ip should be included in the PDF, while one with 
PgT > [if spoils the cut fl3]). However, this is not true at NLO. DGLAP evolution is written in 
terms of parton virtualities q 2 = q T /(l — z), where z is the fraction of parent parton momentum 
carried by the next (in this case, final) parton. So, a relatively soft gluon with p gT < [i F may 
correspond to q 2 > [i 2 F , and thus be assigned to the matrix element. However, this will happen 
mainly for large z close to 1, that is, in a situation where the emission of an additional (and 
now soft) gluon does not change the mass mhard of the 'hard block' too much; and thus does 
not smear out the low x of the gluon sampled by the process. 



4 Procedure to calculate the effect of the proposed cuts 

To demonstrate how effective the proposed cuts will be in determining the gluon PDF at 
small x, we calculate the expected NLO cross section using two different recent sets of parton 
distributions, namely MSTW08 [2] and CT10 [3J, for an LHC energy of y/s = 7 TeV. To 
implement the cuts we adapt the subroutines for the matrix elements squared of LO and 
NLO bb production that are given in the public MCFM program [7]. These subroutines use 
the expressions given in [8] for the NLO loop corrections to the 2—7-2 subprocess, and the 
expressions given in [9] for the 2 — > 3 subprocess. Note that the infrared divergences arising 

5 Even for the B — > D/iv decay, exploited in [I], one can restore the full 4- momentum of the f?-meson based 
on three constraints: we know the B- meson mass and its direction (two angles: <fi and 8), since the position of 
the B-meson decay vertex is observed in the detector. 

6 The f?-meson carries more than 80% [6] of the original b quark momentum. In this paper, we present 
results for bb production. 
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from the emission of very soft gluons are regularized by the so-called 'plus' prescription, which 
means that the singularities in the integrands, as the momentum fraction z — > 0, are tamed by 

I dzf(z) 
Jo 

The 'plus' prescription is well justified at very low z, and corresponds to the Bloch-Nordsieck 
procedure for soft gluon radiation and for the inclusive cross section, where the f(0) terms, 
added in the virtual loop and in the real NLO contributions, cancel each other. However, in 
the 2 — > 3 matrix element not only soft gluons are emitted. If z is not small, the cancellation 
may be spoiled by our cuts. Therefore, to calculate the real 2—7-3 contribution we impose 
the restriction that the virtuality of any external (gluon or 6-quark) line, after gluon emission, 
must be larger than the factorization scale fip. All contributions with smaller virtualities are 
included either in the incoming parton distributions or in the quark fragmentation functions. 
The remaining 2—7-3 contribution, with exactly the same kinematics, due to the additional 
f(0) term, goes to the NLO loop correction to cancel the corresponding (unphysical, if z is 
not too small) term in the loop correction to the 2—7-2 subprocesd3 coming from the 'plus 
prescription'. 

As default parameters we take the renormalization scale to be fiR = mi, = 4.75 GeV, and 
the factorization scale to be [ip = 2 GeV. Also we take k = 2 GeV in ([3]) for the cut on the 
vector sum of the transverse momentum of the outgoing b quarks. 
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5 Factorization scale dependence 

To illustrate the dependence of the predictions for bb production on the choice of the factor- 
ization scale, fiF, we evaluate the cross section for the production of b and b quarks with both 
of their pseudo-rapidities in the interval 5 < r/x,2 < 6, first using \xp = 2 GeV, and then for 
[Ap — 4 GeV. We repeat the exercise for the interval 2 < rji^ < 3. We use the CT10 NLO set of 
partons [3], which are available for Q > Qq = 1.3 GeV. For both choices of rapidity intervals, 
the cross section calculated with the higher scale, \ip = 4 GeV, is about 3-4 times larger than 
that calculated with /ip = 2 GeV. Such a strong factorization scale dependence is due to the 
behaviour of the incoming parton densities. In the small x domain, relevant for the LHC, the 
summation of the double logarithmic terms, 

E nCn (a s ln(l/x)ln(np/Q 2 Q )) n , (6) 

in the DGLAP evolution, leads to an 

exp U (AN c a s /Tx) ln(l/z) ln(/4/Q§)) (7) 

7 The contribution corresponding to /(0) was calculated taking the matrix element for very soft gluon emis- 
sion, keeping the momenta of incoming partons and the outgoing 6-quark fixed. After this, the 1/z singular 
factor was replaced by the corresponding 1/z factor for the 2 — > 3 event which satisfies all of our proposed cuts. 
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growth of the gluon density with increasing /j,p. The exponential growth comes from the 
sum over the possibilities of emitting different numbers of gluons. The growth cannot be 
compensated by the 'hard' matrix element, which at NLO level, allows for the emission of only 
one gluon. This double-logarithmic effect is the main source of the strong factorization scale 
dependence of the predictions for the single 6-quark inclusive cross section. 

On the other hand, if we choose a 'large' value of the scale, \ip > k Q , then we invalidate our 
proposed 'pr' cut (j3J). Recall that an integrated parton density at a scale fip includes the effects 
of all partons with transverse momenta k t < jMp] and the transverse momentum of an incoming 
parton with a 'large' k t will spoil the px balance in Q. To control the transverse momenta 
of the incoming partons we may use unintegrated parton distribution^], f g (x, k t , Hp), an d then 
integrate them over all kt < ko- These distributions can be obtained to NLO accuracy from 
the conventional integrated PDFs following the prescription of Ref. [10]. However, recall that 
integrated PDFs are not available at low scales, k, less than Q . Therefore we replace that part 
of the integral over the unintegrated PDFs with k 2 = k 2 /(l — z) < Qq by the 'integrated' value 
xg(x,Qo)T(Qo, fip), where the Sudakov factor, T, accounts for the probability not to emit an 
extra parton, and thus not to enlarge k t , during the evolution from Qo to \x F . The T factor is 
given by 

T(Q ,M = exp f- ^L^fl f dz <o+nlo ( A (g) 

where the precise form of the splitting function P^® +NLO (z) is given in [TU] . 

With such a procedure, for fixed ko, the main double-logarithmic effects are correctly in- 
cluded via the unintegrated PDF, f g , while the single-log dependence of f g on fi F is mainly 
compensated by the NLO loop corrections in the 'hard subprocess' cross section. The net effect 
of this procedure, based on 'unintegrated PDFs', is a great reduction in the dependence on the 
choice of the factorization scale. For example, changing the scale fip from 2 to 4 GeV now leads 
to less than 20(35)% decrease in the prediction of the bb cross section in the intervals 2 < rj < 3 
(and 5 < i] < 6), rather than the factor of 4 (3) increase, see Fig. [4] which we will introduce 
below. Note that within a 30% accuracy, the result obtained using the unintegrated PDFs 
coincides with that calculated in the conventional NLO collinear framework with fip = ko- 

6 Renormalization scale dependence 

The dependence of the cross section on the choice of renormalization scale, arises from 
the running QCD coupling, a s (fiu). In general, the variation of the QCD coupling should be 
compensated by the logarithmic terms, 1ii(/xr//xf), in the NLO (and higher order) virtual loop 
corrections [IT] . Unfortunately, in the region of interest, say, ftp — 2 GeV and fiR ~ (1 — 2)ra&, 
the corresponding NLO contribution, that is the factor ln(/iR//xp), is practically constant, 

8 We also included the contributions from the incoming quarks, which are, however, negligibly small in the 
low x region of interest. 
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and does not compensate the variation of the LO term, |.M LO | 2 oc o: 2 (//r). Therefore the result 
calculated with = 2m& turns out to be about 30% smaller than that with /i/j = m^. Taken 
together with the uncertainty in the value of the mass, m^, of the 6-quark, this leads to an 
unavoidable uncertainty (~ 50%), at NLO level, in the normalization of the bb cross section^. 

Recall that at these scales the NLO corrections are rather large, larger than the size of the 
LO contribution. Thus it is not evident that the much more complicated NNLO calculation 
will improve the accuracy significantly 



7 Results for bb production after cuts 

Nevertheless, in spite of the uncertainties in the normalization that we discussed above, the 
expected ratio of the cross sections measured in different rapidity intervals is quite stable and 
is driven entirely by the x behaviour of the gluon. We illustrate below how this enables LHC 
bb data to determine the shape of the gluon PDF in the interval 10~ 5 < x 10 -3 . 

Recall that after imposing the cuts of (E} and (@J, the variation of the mass, mh ar d, created 
in the hard subprocess, is strongly limited. The contribution of the 2 — > 3 subprocess never 
exceeds 40% of the whole cross section; typically it only amounts to about 1/3. Moreover, this 
2 — > 3 contribution arises from relatively soft gluon emission, which does not change m hard very 
much. An important consequence is that these bb events in different intervals of pseudo-rapidity, 
1], sample the gluon in rather narrow intervals of x, which allows a precise study of the shape 
of gluon x distribution. 

We illustrate this in Fig. |3j It shows the distributions of the momentum fraction X2 (where 
x 2 < Xi) carried by the gluon in bb events if cuts ()3]) and (jl]) are imposed. The predictions 
are shown for four different intervals of the pseudo-rapidity, rj = 2.5, 3.5, 4.5, 5.5 in flU), 
corresponding to the b and b quarks both having rapidities in the intervals (2,3), (3,4), (4,5) 
and (5,6) respectively. In each case we show the results obtained using MSTW08 and CT10 
parton sets. We take \if = 2 GeV and k Q = 2 GeV. We see, for instance, that the bb events, 
selected by the rapidity cut (3,4), sample the gluon in quite a narrow range about x = 10 -4 . 

Fig. [3] is obtained using integrated conventional PDFs. Fig. H] shows the effect of using the 
procedure based on unintegrated PDFs (as introduced in Section [5]). The figure, obtained via 
CT10 partons, illustrates two effects. First, it shows the reduction in sensitivity to a change 
in factorization scale (from fj,p — 2 to 4 GeV) if unintegrated partons are used. Second, it 
shows that for our default choice, \if = 2 GeV, the predictions obtained using integrated and 
unintegrated PDFs are reasonably similar. 

We have emphasized that, despite normalization uncertainties, the expected ratio of the 
cross sections measured in, say, the pseudo-rapidity intervals (2,3) and (5,6) is quite stable and 

9 For our computations, we choose the conventional value, \xji = mi, = 4.75 GeV [Bl [T^] which was used to 
describe bb production at the Tevatron. 
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Figure 3: The bb cross section as a function of the momentum fraction of the proton carried by the 
slowest gluon, after cuts (T3]) and fl4j) have been imposed, for four different pseoudo-rapidity intervals 
taking r\ = 2.5, 3.5, 4.5, 5.5 in fl4j). In each case we show the predictions obtained using the NLO 
MSTW08 [2] and CT10 [3] parton sets. 

is driven entirely by the x behaviour of the gluon. If the gluon had a pure power behavior, 
xg(x,fip) oc x~ x , then we would observe a flat rj dependence. This would follow since 



A non-trivial, that is non-flat, 77 behaviour of the cross section will reflect the curvature (or a 
deviation from the power law) of the x dependence of the gluon, and may be used to distinguish 
between the different sets of 'global' PDFs. 

For example, in the case of CT10 integrated/conventional partons the expected ratio for 
Hf — 2 GeV is 2.86 (± 15% if we use the unintegrated PDF with fip = 2 or 4 GeV), while for 
the case of integrated MSTW08 partons the analogous ratio is 3.95. 

Therefore a study of the 77 dependence of the bb events, selected by the cuts, is a valuable 
way to study the small x behaviour of the gluon. In Fig. |S(a) we present the rj dependence of 
the cross section expected for six different NLO sets of PDFs, to demonstrate the sensitivity of 
such a method to the small x behaviour of gluons. For this plot, X2 varies from X2 ~ 3 x 10~ 4 
(corresponding to rj = 2.5) to 12 ~ 2 x 10 -5 (corresponding to 77 = 5.5). This corresponds to 
the pseudo-rapidity range relevant to the LHCb experiment. In this small x region the PDFs 
are unconstrained by existing data. In particular, we note that the gluon distribution is well 
determined at Q 2 = 5 GeV 2 to within 20% (typically less than 10%) only if x is in the range 
10~ 3 < x < 0.3, see Fig. 16 of [2J. Since the gluons are unconstrained at very low x, we should 
anticipate the large spread of the predictions, shown in Fig. EJ^a), for the rj dependence of the 
cross section after the cuts are imposed. Indeed, if the errors on the PDFs were to be included, 
then the various predictions would overlap. 
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Figure 4: The bb cross sections as a function of the momentum fraction carried by the slowest 
gluon, after cuts (El) and (0} have been imposed, predicted using integrated and unintegrated PDFs 
for two choices of factorization scale, fip — 2 and 4 GeV. The prefix "u" indicates unintegrated 
PDFs are used. The left and right plots correspond to the rapidity intervals (5,6) and (2,3) specified 
by taking r? = 5.5 and r] = 2.5 in fllj. The \if = 4 GeV prediction obtained from integrated partons 
is not shown for the latter interval since it about 4 times higher than that for jj F = 2 GeV. CT10 
NLO partons [3] are used. 



Recall that there is an overall normalization scale uncertainty of about 50% (see Section 
EJ). Even allowing for this, the observed bb cross section, da/drj, after cuts, will be provide 
valuable information about the gluon in this completely unexplored low x domain. However, 
the normalization uncertainty does not affect the shape of the prediction of the cross section 
versus rj, so the bb data are capable of yielding even more precise information. Simply for 
illustration, we compare in Fig. MJo) the shapes of the predictions of all the parton sets after 
they have been normalized to 1 at 77 = 2.5. Of course, the present huge PDF errors mean such 
a comparison cannot distiguish between them. Rather the bb data will determine the gluon 
distribution to 50% and its shape to much better accuracy. Clearly the shape has increasing 
discriminatory power as 77 increases. 

We emphasize again that we have used the integrated PDFs obtained from the various 
'global' analyses at face value, simply to illustrate that bb data at the high values of 77, accessible 
at LHCb, offer a powerful probe of the gluon distribution at small x, provided that the cuts given 
in ([3]) and (jl]) are applied. The LHC data will be able to probe a low x, low Q 2 domain well 
beyond the range of the data fitted in the 'global' PDF analyses. Of course, in this domain, 
we have no reason to trust any of the predictions obtained from 'extrapolated' PDFs. We 
merely show the curves obtained from six different PDF sets to illustrate the potential ability 
to constrain the gluon in the 10 -5 < x < 10~ 3 domain by measuring bb events at high 77, with 
cuts (j3J) and @ imposed. 
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Figure 5: (a) The r\ dependence of the bb cross section, after the cuts (TJD and (0} have been 
imposed, obtained using six different NLO sets of partons MSTW08, CT10, NNPDF21, ABKM09- 
4, HERAPDF01 and GJR08VF (b) the predictions normalized to one at rj = 2.5. 



8 Conclusion 

At LHC energies, the bb cross section predicted within the conventional NLO collinear approach 
has huge uncertainties (up to a factor of 4) arising from the variation of the factorization 
scale, \ip. This uncertainty can be strongly reduced by selecting events where the transverse 
momenta of the two 5-mesons balance each other to some accuracy. If the sum of the momenta 
\p 1T + p 2 r\ < ^0) then the scale \ip ~ k . This offers the possibility to use bb production at the 
LHC, particularly in the LHCb experiment, to study the x-dependence of the gluon distribution 
down to x ~ 1CT 5 at rather low scales, /i 2 ~ 4 GeV 2 , where the present HERA and Tevatron 
data do not constrain the behaviour of the parton densities, and where different parton analyses 
propose quite different gluons. 

We have considered the renormalization and factorization scale dependences of the bb cross 
section after the cuts ([3]) and (j3J), on the transverse momenta and rapidities of the 5-mesons, 
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have been imposed. In this way, we have demonstrated how such events may determine the be- 
haviour of the gluon distribution down to x ~ 10~ 5 . In Fig. [5] we showed the differences obtained 
using six different set of PDFs^l simply for illustration, bearing in mind that extrapolations of 
existing PDFs into this domain are unreliable. 
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